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Abstract. In 2010, several experimental groups obtained compelling evidence for spin-triplet 
supercurrent in Josephson junctions containing strong ferromagnetic materials. Our own best 
results were obtained from large-area junctions containing a thick central Co/Ru/Co "synthetic 
antiferromagnet" and two thin outer layers made of Ni or PdNi alloy. Because the ferromagnetic 
layers in our samples are multi-domain, one would expect the sign of the local current-phase 
relation inside the junctions to vary randomly as a function of lateral position. Here we report 
measurements of the area dependence of the critical current in several samples, where we find 
some evidence for those random sign variations. When the samples are magnetized, however, 
the critical current becomes clearly proportional to the area, indicating that the current-phase 
relation has the same sign across the entire area of the junctions. 



The interplay between magnetism and superconductivity has been a fertile ground for new 
physics over several decades. On the one hand, bulk materials exhibiting both forms of order 
are rare, although there is recent progress in discovering new materials with exotic properties 
[UE1E]- On the other hand, combining conventional superconducting (S) and ferromagnetic (F) 
materials relaxes the constraints, leading to a wide variety of possibilities. The interplay between 
the two is governed either by the magnetic field [I] or the exchange field [5j; in this work we 
restrict our attention to the latter. The essential physics arises from a modification of Andreev 
reflection - the microscopic process by which electrons enter or leave S. A conceptually simple 
picture is to consider two electrons belonging to a Cooper pair, which happen to cross the S/F 
interface. Because they enter opposite spin bands in F, they have different Fermi wavevectors, 
hence the pair acquires a center-of-mass momentum perpendicular to the S/F interface [6]. 
That, in turn, leads to oscillations in the pair correlation function, which manifest themselves as 
oscillations in the T c of S /F bilayers [7] , inversion in the proximity- induced density of states on 
the F side of S/F bilayers [8], and oscillations between "0" and 'V states in S/F/S Josephson 
junctions [10] . The oscillatory phenomena associated with S/F hybrids are fascinating and 
have potential for some applications; on the other hand, they are very short-range [5], decaying 
on a length scale governed by the exchange energy or by the mean free path in F. 

A new chapter in S/F physics was opened in 2001 with the prediction that spin- 
triplet pair correlations can be induced in S/F hybrids made from conventional spin-singlet 
superconductors, in the presence of specific types of magnetic inhomogeneity that include non- 
collinear magnetizations [12| \T3[ H3] • Since the two electrons from a spin-triplet pair with 
m = ±1 (where m is the component of spin angular momentum along the magnetization direction 
of F) live in the same spin band in F, they are not subject to the exchange field. As a result, 



such spin-triplet pair correlations are long-range in F, decaying on a length scale governed by 
the temperature or by spin-flip or spin-orbit scattering. The experimental search for spin-triplet 
pairs produced a few promising hints a few years ago |15|. [T6] , but conclusive evidence for their 
existence came only last year when several groups reported long-range supercurrents in S/F/S 
Josephson junctions [TTJ, [18] US EQ]. Our own results [T7] were based on junctions of the form 
S/N/F'/N/SAF/N/F"/N/S (see Fig. 1), where SAF is a Co/Ru/Co "synthetic antiferromagnet" 
|21j . N is Cu, and F' and F" are thin ferromagnetic layers (Ni in this work). Conversion of spin- 
singlet to spin-triplet pairs in our samples can be viewed as taking place in two steps [Hj: 1) 
The two electrons of a spin-singlet Cooper pair leaving the bottom S layer acquire different 
phase factors in the F' layer, which generates the m = triplet component; 2) Rotation of the 
m = triplet component into the basis defined by the SAF magnetization generates the m = 1 
and m = — 1 triplet components in the new basis, as long as it is non-collinear with the F' 
basis. Conversion back to spin-singlets occurs between the SAF and F" layers as the inverse 
process. This type of junction, with three ferromagnetic layers, had been previously proposed 
theoretically by Houzet and Buzdin [22], but with a simple F layer in place of the SAF, and 
without the crucial inner two N layers which serve to magnetically decouple the F' and F" layers 
from the Co layers. The critical supercurrent in these samples hardly decays as the total Co 
thickness varies from 12 to 30 nm |17| . in sharp contrast to similar samples without F' layers 
where the supercurrent decays very rapidly with increasing Co thickness [23] . 
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Figure 1. Schematic diagram of the 
Josephson junctions used in the this work, 
shown in cross-section. The purpose of 
the Au layer is to prevent oxidation during 
processing; at low temperature it becomes 
fully superconducting by the proximity effect. 



The results described above [17] immediately raise several questions, only some of which 
have been addressed by subsequent work [23]. One of the most important outstanding questions 
concerns the current-phase relation of the junctions, or more precisely, the current-phase relation 
as a function of position across the lateral dimensions of the junction. At issue is the theoretical 
prediction that S/F'/F/F"/S or S/F'/SAF/F"/S junctions (we omit explicit reference to the N 
layers for simplicity) can be either "0" or 'V junctions depending on the relative orientations of 
the magnetizations in the three ferromagnetic layers |22 [ I25 [ [26 [ I27|. (In the weak-coupling limit, 
a 0-junction has the current-phase relation I s = I c sin(</>), whereas a 7r-junction has the current- 
phase relation I s = I c sin(0 + ir) = —I c sin(0), where (ft is the gauge-invariant superconducting 
phase difference between the two S electrodes.) Looking through the junction in the direction 
of current flow, the three magnetizations can either rotate in the same direction from one to 
the next (with either chirality), or the rotation from the F' layer to F can be in the opposite 
sense with respect to the rotation from F to the F" layer. (We define the direction of rotation 
using the angle that is less than ir.) In the former case, the junction is a 0-junction, whereas in 
the latter case it is a 7r-j unction [27]. Since our samples have relatively large lateral dimensions 
(10, 20, or 40 /im in our published work [TT] ) , and since a typical Co domain size is only on the 



order of 3/xm in our sputtered films in the relevant thickness range [28], we expect our samples 
to contain many domains which produce a random pattern of and n junctions. In that case 
one would expect the typical critical current of a sample to be proportional to the square root 
of the number of domains, hence proportional to the square root of the junction area. This is 
quite different from the usual case where the current-phase relation is uniform across the entire 
area of the junction, and the critical current is proportional to the junction area. The main goal 
of this paper is to determine experimentally which type of area scaling occurs in our samples. 

Our sample fabrication procedure is described in detail in our previous publications 
[29\ 117] . All of the samples discussed in the present work contain two 1.5-nm-thick Ni 
layers as the F' and F" layers as shown in Figure 1. Photolithography and ion milling are 
used to define Josephson junctions with circular cross-section. Most of our published data were 
obtained from junctions with diameters of 10 or 20 /im. Taking into account sample-to-sample 
fluctuations, we were not able to determine the area scaling law from our previous data sets 
\17\ [23] . Since the goal of the present work is to measure how the critical current scales with the 
junction area, we fabricated junctions with diameters of 3, 6, 12, and 24 /im on each substrate. 
Although the success rate of the largest junctions is low, we obtained data from 3, 6, and 12 /im 
diameter junctions on several different substrates. All of the measured junctions display I — V 
characteristics obeying the standard overdamped form: V = Rn{I 2 — I 2 ) 1 ^ 2 , where Rn is the 
normal state resistance. 
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Figure 2. Josephson junction 
critical current vs. magnetic 
field applied in the plane of 
the substrate - perpendicular 
to the current direction. In 
the virgin state (panels a-c), 
the resulting Fraunhofer pat- 
terns vary in quality. After 
the same samples are magne- 
tized (panels d-f), the quality 
of the Fraunhofer patterns im- 
proves and the maximum crit- 
ical currents are much larger. 
The data shown were taken at 
T = 4.2K from junctions fab- 
ricated on a single substrate, 
with junction diameters of 3, 
6, and 12 /im for the upper, 
middle, and lower panels, re- 
spectively. 



For every Josephson junction we fabricate, we measure the dependence of I c with magnetic 
field applied perpendicular to the current direction, i.e. in the plane of the substrate. The left 
panels of Figure 2 show the resulting "Fraunhofer patterns" of a typical set of junctions in the 
virgin state, fabricated on a single substrate, with diameters of 3, 6, and 12 /im. The quality of 
the virgin-state Fraunhofer patterns varies, which is not surprising given the amount of magnetic 
material contained in the junctions [30} 129]. If the SAF is working effectively, then the Co layers 
produce zero net magnetic flux in the junctions, because a given domain in one Co layer is 
antiparallel to that in the other Co layer [23]. Deviations from ideal Fraunhofer patterns are 
most probably due to the two Ni layers. Figure 2 also shows Fraunhofer patterns of the same 



three junctions after they have been subjected to a large in-plane magnetic field of 1800 Oe. (The 
large field is then removed to obtain the low-field Fraunhofer pattern.) One immediately notices 
three features in the data: 1) the quality of the Fraunhofer patterns is improved; 2) the central 
peak of the pattern is shifted to a small negative field; and 3) the maximum critical current 
determined from the largest peak in the pattern is much larger than in the virgin state data. All 
three features can be explained by two hypotheses: i) the Ni layers become fully magnetized in 
the direction of the applied field; and ii) the SAF undergoes a "spin-flop" transition, whereby 
the two Co layers end up perpendicular to the direction of the applied field |3T1 [32]. The net 
result is that the Ni and Co magnetizations end up perpendicular to each other, which optimizes 
generation of the spin-triplet pairs \Z2 \ \2<5 \ |2"T]. Meanwhile, the central peak of the Fraunhofer 
pattern is shifted to a small negative field where the magnetic flux due to the external field 
exactly cancels the magnetic flux due to the Ni magnetization [333 122] • Further evidence for the 
occurrence of the spin-flop transition in our samples will be provided elsewhere |34] . 



Figure 3. Critical current times 
normal-state resistance (I c Rn) as a 
function of magnetizing field H app for 
the three different diameter Josephson 
junctions shown in Figure 2. The large 
increase in I c Rn occurs in the vicinity 
of the coercive field of the outer Ni 
layers, which is approximately 800 Oe 
for the Ni thickness of 1.5nm in these 
samples. All data were taken at T = 
4.2K. 

The evolution of the I c enhancement as the samples are magnetized is shown in Figure 3. 
The figure shows I c Rn as a function of H app , the magnitude of the applied field, as that field 
is increased in discrete steps. After each step, the full Fraunhofer pattern is measured near 
zero field; the quantity plotted is from the value of I c obtained from the peak of the resulting 
Fraunhofer pattern. For all three samples, I c starts to increase when H app exceeds about 800 
Oe, and attains a plateau when H app reaches about 1800 Oe. Data from samples with different 
Ni thicknesses [M] show that the values of H app where I c increases vary with Ni thickness in a 
manner consistent with the coercive field of thin Ni films (larger coercive field for thinner films). 
The field at which the spin-flop transition occurs cannot be determined from these data; all we 
know is that it is less than 500 Oe - the coercive field of our thickest Ni layers measured. 

Returning to the question of how I c scales with junction area, Figure 4 shows a plot 
summarizing data from 18 Josephson junctions on 6 separate substrates. We plot I c Rn in both 
the virgin state and after the samples are fully magnetized, as a function of junction diameter. 
Several systematic features of the data stand out. First, the data obtained from the magnetized 
samples are remarkably consistent over all junction sizes and all substrates. Since Rn is inversely 
proportional to the junction area A, the fact that I c Rn is independent of area implies that I c 
scales linearly with junction area in the magnetized state, i.e. I c oc A. This means that the 
supercurrent is adding coherently across the area of the junction. That is expected when the 
two Ni layers are magnetized in the same direction across the whole junction area. For example, 




if we define the direction of H app as 9 = 0, so that the two Ni magnetizations point along 9 = 0, 
then the lower and upper Co layers should point either along 9 = ir/2 and —n/2, respectively, 
or the converse, due to the spin-flop. In either case, the sense of rotation from the bottom Ni 
layer to the lower Co layer will be the same as the sense of rotation from the upper Co layer to 
the top Ni layer. This would be true even if the Co magnetizations chose an angle different from 
±7r/2, as long as that angle is not and the two Co magnetizations point in opposite directions. 

The virgin state data shown in Figure 4 show more variation from substrate to substrate. 
Before discussing those data, we first mention some variations in the ion milling process we used 
to fabricate the samples. To produce Josephson junctions with a well-defined lateral area, it is 
only necessary to ion mill through the upper 25-nm-thick Nb layer (see Fig. 1). We generally 
continue milling through the two Cu layers, the F' layer, and part-way through the upper Co 
layer to ensure that there is no Nb left outside of the junction area. Of the samples reported 
in this study, samples #9, #12 and #14 were ion milled part-way through the upper Co layer. 
In these samples the lower Co layer is unaltered by the ion milling, and should maintain the 
domain structure of the continuous film. For these samples, panel (a) of Figure 4 shows that 
the virgin-state I c Rn decreases with increasing junction diameter D, with an almost inverse 
linear dependence (except for the 12^m pillar of sample #12). That means that I c Rn varies 
nearly inversely with the square root of the junction area, I c Rn 3c A -1 / 2 , which means that I c 
is proportional to the square root of the area, I c oc A 1 / 2 . That is what one would expect if the 
junctions consist of multiple domains contributing a random distribution of and tt junctions. 
In this scenario, the large increase in I c that occurs when the samples are magnetized is largely 
due to a transition from I c oc A 1 / 2 scaling to I c oc A scaling. Optimization of the angle between 
the Ni and Co magnetizations to tt/2 also plays a factor [34| : we cannot determine from our 
data which factor plays the larger role. 




Figure 4. Summary of I c Rn vs. Josephson junction lateral diameter for all samples measured, 
in both the virgin state (solid symbols) and after magnetization (open symbols). Each different 
symbol shape and color refer to a different substrate. 

Panel (b) in Figure 4 shows data from several substrates that were ion milled further - 
either part-way through the lower Co layer (samples #15 and #16) or in one case all the way 
through the lower Co layer (sample #7). In those samples the scaling of the virgin state data 
is less clear, and may even be consistent with I c oc A scaling. It is possible that further ion 
milling of the central Co/Ru/Co SAF and even the bottom Ni layer (in the case of sample #7) 
imposes constraints on the domain structure of the ferromagnetic layers, so that the domain 
magnetizations can no longer be described as pointing in random directions. 

The large sample-to-sample fluctuations in I c Rn seen in the virgin state data in Figure 4 



are strongly correlated with the quality of the Fraunhofer patterns. For example, the 12 /xm 
junction of sample #12, which clearly deviates from the I c oc A 1 ! 2 scaling proposed above, had 
a particularly nice Fraunhofer pattern - much nicer than those from the 3 and 6 //m junctions 
on the same substrate. Similar correlations were observed in other samples: e.g. the 6 //m 
junction of sample #16 had a nicer Fraunhofer pattern than did the 3 or 12 /im junctions on 
that substrate. 

In conclusion, we have observed strong I c oc A scaling in S/F'/SAF/F"/S Josephson junctions 
after they are magnetized, demonstrating coherent superposition of supercurrent over the whole 
junction area. In the virgin state, there is evidence for I c oc A 1 / 2 scaling in some samples, but 
not in others. Large variations in the quality of the Frauhofer patterns acquired from the virgin- 
state samples partially mask the area scaling. We plan to repeat this study using Pdi-^NLj 
alloy in place of Ni for the F' and F" layers, since our previous work showed that PdNi layers 
do not distort the virgin-state Fraunhofer patterns [29]. 
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